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1 Introduction 

In this paper we consider the following nonlocal problem with integral con- 
dition: 

^ + Am = 0, t e [0, T] 

rT (1) 

u(0) + / w{s)u{s)ds = uq, 



where w{s) > is a given function, uq G X. The operator A with the 
domain -0(^4) in a Banach space X is assumed to be densely defined strongly 
positive (sectorial) operator, i.e. its spectrum S(A) lies in a sector of the right 
half-plane with the vertex at the origin, while its resolvent decays inversely 
proportional to \z\ at the infinity (see estimate ([6]) below). 

Inhomogeneous problem related to ([1]) can be reduced to homogeneous 
one by change of function in the following way. If we have 

^ + Av = f{t), te[o,r] 

rT (2) 

v{0) + / w{s)v{s)ds = uo, 
Jo 

with /(t)- vector- valued function in the Banach space X then putting v(t) = 

u(t) + Vi(t), where 

t 



we obtain the following problem for u(t) 

— + Au = 0, t e [0, T] 

(JjO 

r 

m(0) + / w{s)u{s)ds = Mo — $, 
Jo 

where 

T 

$ = / w{s)vi{s)ds. 





Note that exponentially convergent numerical approximation for vi (t) was 
developed in |1], [3]. So, one can use this approximation to obtain vi{t) and 
then to find $. 



It should be noted that recently there were developed various exponen- 
tially convergent methods for problems with unbounded coefficients in Ba- 
nach space [2], [3], |6], [7], [9], [TT]. These problems can be considered as 
metamodels of classical problems for partially differential equations such as 
parabolic elliptic and hyperbolic. 

The aim of this paper is to construct an exponentially convergent approx- 
imation of a solution to problem ([T]). The paper is organized as follows. In 
Section[2]we discuss the existence and uniqueness of the solution as well as its 
representation through input data. A numerical algorithm for the problem 
(^ is proposed and justified in section |3l The main result of this section is 
theorem [1] about the exponential convergence rate of the proposed discretiza- 
tion. The next section H] is devoted to numerical examples which confirms 
theoretical results from the previous section. 

2 Existence and representation of the solu- 
tion 



The solution of ([T]) can be formally represented as follows 

u{t) = e-^'u{0). 
From the integral condition in ([T]) and ([3]) we obtain 

m(0) + / w{s)u{s)ds = Mo, 



(3) 



Therefore, in the case when 
tions see below) we have 







u{0)+ I M;(s)e"^'rfsM(0) =Mo, 



I + Jq w{s)e~ ds exists (sufficient condi- 



So, 



u{0) 



u{t) = e" 



As 



1+ I w;(s)e"^^(is 





Uq. 



-At 



T n-1 



Uq. 



(4) 



Let the operator A in ([T]) be a densely defined strongly positive (sectorial) 
operator in a Banach space X with the domain D{A), i.e. its spectrum S(y4) 
is situated in a sector S 

TT 



S = {^ = po + re'' : r G [0, oo), po > |^| < (^ < ^} . 



(5) 



Additionally, the following estimate for the resolvent is valid 

\\Ra{z)\\ = \\{zI-A)-'\\<-^ (6) 

~r I ■^ I 

outside the sector and on its boundary F^. The numbers po? V^ are called the 
spectral characteristics of A. 

We call a spectral hyperbola the curve Fq : 

To = {-2^(0 = Po cosh C — ifeosinhC : C ^ (~oo, oo), 60 = Po tany?}. (7) 

It has a vertex at (po, 0) and asymptotes that are parallel to the rays of the 
spectral angle S. 

A convenient representation of operator functions is the one through the 
Dunford-Cauchy integral (see e.g. [HE]) where the integration path plays an 
important role. Using the Dunford-Cauchy integral representation and (jl]) 
the solution to problem ([T]) can be represented as 



uit) = / e 



-zt 



2711 ./p 



w{s)e ^''ds RA{z)uodz 
F{z, A)RA{z)uodz, 



I 

if F{z,A) is analytic function inside the integration hyperbola Tj that en- 
velopes Fq. To obtain uniformly convergent and numerically stable algorithm 
we shall modify this integral by changing the resolvent Ra{z) to R\{z) that 
doesn't change the value of integral when uq G D^A""), a > (for the details 
see [3]). 

R\{z) = {zl - A)-' ~ ^-. 

Therefore, one can obtain the following representation for the solution of the 
problem ([1]): 

u{t) = 7^ [ F{z,A)R\{z)uodz. (8) 

We choose the following hyperbola 

F/ = {z{() = aj cosh C — ibj sinh ( : ( E (—00, 00)}, (9) 

for an integration contour that envelopes the spectrum of A, where the values 
of a/, bj are to be defined later. Using this hyperbola, we obtain from ([8]) 

1 /'oo poo 

u{t) = — / F{z{0, A)R\{C)z'iC)uodC = / J-(t, C)dC (10) 

4 



with 

z'{C.) = aj sinh ( — ibj cosh (. 

The next step toward a numerical algorithm is an approximation of ( jTOj) 
by the efficient quadrature formula. For this purpose we need to estimate 
the width of a strip around the real axis where the integrand in (TTOj) permits 
analytical extension (with respect to (). The integration hyperbola Tj will 
be translated into the parametric set of hyperbolas with respect to u after 
changing ( to ( + iu 

T{iy) = {z{(, u) = aj cosh ((^ + iu) — ibj sinh {C, + iu) : ^ G (—00, 00)} 
= {z(C, i^) = a(i^) cosh C — ibi^u) sinh ( : ( E (— C)0, 00)}, 

with 



a(z/) = aj cos i/ + 6/ sin z/ = yaj + bj sin (z/ + 0/2), 

6(i^) = bj cos 1/ — a/ sin 1/ = Wa| + bj cos (i^ + 0/2), 

(j) bj (j) aj 

cos — = — , , sin — 



2 VaJ+6f' 2 v^olTfef 
The analyticity of the integrand in the strip 

Dd, = {(C, ^)-(^ (-00, cx)), |i^| < di/2}, 

with some rfi could be violated if the resolvent or the part related to the 
nonlocal condition become unbounded. To avoid this we have to choose di in 
a way such that for u e [—di/2, di/2) the hyperbola r(z/) remains in the right 
half-plane of the complex plane. For u = —di/2 the corresponding hyperbola 
is going through the point (pi,0), for some < pi < po- For u = di/2 it 
coincides with the spectral hyperbola and therefore for all z/ G {—di/2, di/2) 
the set T{v) does not intersect the spectral sector. This fact justifies the 
choice the hyperbola r(0) = F/ as the integration path. 

Such requirements for F(z/) imply the following system of equations 

ai cos {di/2) + bi sin {di/2) = po, 

bi cos {di/2) — aj sin {di/2) = 60 = Po tan cp, 

aj cos {—di/2) + bj sin {—di/2) = pi, 

it leads us to the next system 

o/ = Po cos {di/2) — bo sin {di/2), 
bi = Po sin {di/2) + bo cos {di/2), 
2aj cos (rfi/2) = Po + Pi- 



Eliminating a/ from the first and the third equations we obtain 

Po cos di — Bq sin di = pi, 
cos((ii + if) 



VPo + ^o' 



I.e. 



di = arccos ( , ^^ ) - ^, (11) 



^0, 



with cosiyC = /" , sin(/) = , 'y . Thus, for a/, 6/ we receive 

di 



«/ = V^o + ^qCos ( y + <^ 



(f + ^)_ ^°^ H°^ (^fe) /' + ^/' 



(12) 



cos 
= Po ^^ = Po 

COS (y9 COS ip 

bi = ypl + blsinl^j + ip 

(d, , \ cos I arccos I /' . ) /2 + (/?/2 
cos(f + (/?) \^ VV^'o+^V 
= Po — = Po • 

cos (f COS !p 

For a/ and hi defined as above the resolvent of the operator A is analytic in 
the strip Dd^ with respect to w = C,+iv for any t > 0. Note, that for pi = 
we have di = 7t/2 — ip as in |1]. 

Taking into account f ll2p we can similarly write equations for a(z/), 6(i^) 
on the whole interval — -s^ < z/ < ^ 



a{i') = aj cos u + bisinu = J p^ + h^ cos I ir + 'Z' I cos(i/) 

+ V Po + ^0 sin ( y + ^ j sin(z/) = y^p^ + 52 ^qs ( y + y? - z/ 
6(z/) = hi cos z/ — a/ sin 1/ = wpg + 6q sin I — + p\ cos(z/) 

^Po + ^0 cos ( y + (^ j sin(z/) = y'p^ + feg gin ( y + ^ - z/ 



Pi < a{v) < Po, 60 < K^) < \ bl + Po - P? 



with di, defined by (ITTl) . 

Now, let us establish condition on w{s), when expression 



1 + / w{s)e~'"'ds 





related to nonlocal condition dose not become zero inside the integration 
hyperbola F/. 



w[s)e 



-^(^)^rf. 



> 1 



wisje 



~<^>ds 



> 



>l-\\wis)\\c[o,T] I e-^"^™^- 
/o 



> 



aj cosh(C) 

_ l|w(g)||c[o,r] 
aj 



> 



Therefore, we have 



w(sie 



-^^<>ds 



< 



1 



\h{s)\\clO.T] 
0.1 



<Ci, 



in the case when 

lk(s)||c[o,T] < a/, 

where a/ is defined in flT^ . 

So, we can summarize all of the above in the following lemma. 



(13) 



Lemma 1 Let A he a densely defined strongly positive operator. If the con- 
dition ( !T3|l is valid then there exists a unique solution to problem ^ that 
can he represented by (jS]). 

More rough estimate than (IT3l) is 



\w{s)\\c[0,T] < 



r 



(14) 



that can be easily obtained from estimate 



W[S]e 



-''^<>ds 



< \\w(s) 



\C[0,T] 



T. 



Further, let us establish conditions for the existence of the solution to ([T]) 
in the case when the operator A is self-adjoint positive definite. To achieve 



that we have to choose di in a way that for v G (— di/2, di/2) the hyperbola 
T{v) remains in the right half-plane of complex plane. For v = —di/2 the 
corresponding hyperbola turns into the line parallel to the imaginary axis. 
For i/ = di/2 it coincides with the ray that lies on the real axis having a 
vertex at po- These requirements imply the following system of equations 

aj cos {di/2) + bj sin {di/2) = pq, 
hj cos {di/2) — aj sin {di/2) = 0, 
aj cos {—di/2) + bj sin {—di/2) = 0, 

which has the solution 

The condition ( 1131) then becomes 

||w(s)||c[o,r] < -|, (15) 

that is sufficient condition of existence solution to ([1]) in the case of self- 
adjoint positive operator A. 



3 Numerical algorithm 

First of all we approximate integral in denominator using exponentially con- 
vergent quadrature. For this aim we use Gauss quadrature rule. So, 

rT ^ rp 

1= w{s)e-'^^^'ds ^ V -UjwiQe-'^'^^^^ = /„, (16) 

e, = |(^, + i), 

where 6j- are a set of n + 1 roots of the Legendre polynomial Pn+i(x) and 
Uj- a are a set of weights related to the Gauss quadrature rule. Note that 6j 
and Uj can be precomputed using fast algorithms (see [TO]). 
Therefore we obtain from fITOl) 

w(t) ^ w„(t) = — / F„(^(C), A)i?],(C)^'(C)%c?C = / J-„(t, C)dC, (17) 



where 



For the error estimate we have 
1 1 



1 + / 1 + /„ 
Due to (1131) or (IHD we have 



In -I 



|1 + /| 
For the second muhipher we obtain 



;i + /)(! + /„ 



<C. 



\l+Ir, 



< 



1 



J2 ujjw{Qe-''^Oi^ 



j=0 



< 



1- 



^ll"'(s)llc[0,T] 



< 



3=0 



< 



1 



< c = const, 



1 -T||w;(s)||c[o,T] 
in the case when (IT4D is vahd. Therefore we have 



<c|J„-J|. 



1 + / 1 + 4 



Exponential function e~^* is analytical on s in all complex plane. So, smooth- 
ness of integrand in I is defined by w{s). Using theorem 19.3 from [TU] we 
have that if w(^{s + 1)) is analytic in [—1, 1], analytically continuable to the 



open Bernstein ellipse where 

|/-/nl < 



144Mp-2" 
35(p2-l)' 



< M then 

n> 2. 



(19) 



If w{s) and its derivatives through w^^ ^' are absolutely continuous and w^^' 
is of bounded variation V then 



|/-/„|< 



32^ 



n > 2i/ + 1. 



157ri/(n - 2z/ - l)2'^+i' 
Supposing uo G D(yl°), < a < 1 it was shown in [3] that 

||e--(C)V(C)i?i(C)Mo|| < (1 + M)K— f — ^ e-Mcoshc-a|Clp" 

C e M, t > 0, 



(20) 



(21) 



where K is a constant that depends on a, M is a constant from resolvent 
estimate (E]). 

The part responsible for the nonlocal condition in ( 1T7|) is estimated by 
(TT8|) . Thus, we obtain the following estimate for J^(t,^): 



C{^,a) 



|^n(t,C)|| < C(v.,«)e-'^^*^°^^'^-"l'^lp"Mo| 

_ {i + M)Kcbi r 2Y 

2iTai \ai 



C G M, t > 0. 



(22) 



We approximate integral flT7|) by the following Sine- quadrature formula 
Pi: 



TV 



Un,N 



{t) = hY, j^^{t,z{kh)), 



(23) 



k=-N 



with the error 



\'r]NiJ^n,h)\\ = \\Unit) -M„,Ar(t)|| 



< 



U, 



,{t)-h Y, Mt,4kh)) 



k=—oo 

I Q-TTd/h 



+ 



h J2 Mt,4kh)) 

\k\>N 



< 



47rsinh(7rd//i)""^""^'(^'*) 



-C{^,a)h\\A''uo\\ Y, 6" 



ait cosh kh—akh 



k=N+l 



Here ii^{Dd) is a space introduced similarly to [8] in [3] of all vector- 
valued functions J-" analytic in the strip Dd ■ Due to [3] 



|e-^(-)V(-)i?]i(>o||Hi(D,) < \\A^uo\\[C4v,a,6) 
+ C+(<^,a,5)] / e-"l«lcie = C(<^,a,5)||A"no|| 



(24) 



with 



C{^, a, 5) = -[C+{^, a, 5) + C.{^, a, 5)], 
C±((^,a,(5) = (l + M)A'tan(^ + ^±- 



2cOS(/9 



° (25) 



4 2 2; \^aocos(f + 1 ±f) 
d = (ii — (5, 



10 



for an arbitrary small positive 6. 

Obviously that in the case of flT^ the part responsible for the nonlocal 
condition is bounded in Dd- Therefore we obtain 

\\MtX)\\mD,)<C{^,a,S)\\A'^uo\\. 
So, we have for the error r]iy{J^n, h) 

WtIaAF h)\\ < ^11^"^°^! /_!l! + -a jt cosh {{N+l)h)-a{N+l)h\ ^^Q) 

WriN^-tn^mS ^ \sinh(7rrf//i)+^ j ^ ^^ 

where the constant c does not depend on /i, A^, t. 
Equalizing the both exponentials gives 



Txdi 



(27) 



this leads us to the following error estimate 

\\TlN{^n.h)\\ < -exp(^-^/nda{N + 1)) \\A"uo\\ (28) 

In the case t > 1 the first summand in the argument of e""^*™'^'^'-^^"''^^^^""''^"'""'^^^ 
in (l26l) contributes mainly to the error. Setting for such case h = ciln N/N 
with some positive constant ci we obtain for a fixed t the following estimate: 

\\VN{J^n,h)\\ < C [e-diN/{c,lnN)^^-c,tajN/2-c,alnN^ P"no||. (29) 

Thus, we have proven the following theorem. 

Theorem 1 Let A be a densely defined strongly positive operator, Uq G 
D[A°'), a G (0,1) and condition flT^ is valid. Then Sine- quadrature fl23|l 
represents an approximate to Un{t). It provides the convergence of exponen- 
tial order uniformly with respect to t > presented by the estimate fl2Sl) for 
the step size h defined in T27}i . The approximation has the convergence rate 
fl29|l for the case t > and h = ci In N/N. 



Remark 1 The integration curve Tj is symmetric with respect to the real 
axis. Therefore z{—kh) = z{kh) and z'{—kh) = —z'{kh). Approximation 
can be rewritten in the form 



Un,N{x) = — - j:„(t, 2(0)) + Re 
Ztii 



sr^ J^n(t,z{kh)) 
^-^ Txi 



which reduce the number of resolvent calculations by factor of two. 

11 



Further we consider the full error estimate 

El = \\u{t) -Un{t)\\ 



I 

2^ 



-z(C)V' 



C'(C) 



/oo 
[j'{t,c)-MtX)]dC 
-oo 

R\{C)uo\\dC 



1 + / 1 + L 



then using OTl) we obtain 



£i 



'l + M)Kbic f 2Y 



< 



27ra/ \ai 

[1 + M)Kbjc [ 2\ 



( — j ||v4"no|| |/ - 41 I e-'^^*'=°^^^-"l^'lciC < 



7ra/ 



||^"^o|||/-/. 
ajj 



-"i^irfC 



= — 11^ Moll |/ - 41 = C \\A Moll |/ - In\ ■ 

For the full error we have 

\\u{t) - Un,Ar(t)|| < £1 + ||?7Ar(J^„, /l) || • 

Therefore we can formulate the main theorem. 

Theorem 2 Let the conditions of theoremUl are valid. Then (125]) represents 
an approximate to u{t). It provides the convergence of exponential order in 
the case when w{t) is analytically continuahle to the Bernstein ellipse. 

4 Numerical examples 

Example 1 Let us consider the problem ([T]) with the operator A defined by 



D{A) = {v{x) e H\0, 1) : v{0) = v{l) = 0}, 
Av = -v"{x) Vt; eD{A), 



(30) 



that generates a homogeneous parabolic equation with boundary conditions 

du{x,t) d'^u{x,t) 



dt dx"^ 

M(0,t) = u(l,t) = 0. 



0, 



Let us supplement this problem with a nonlocal integral condition 



m(x,0)+ / cos{s)u{x,s)ds 
'o 



7r4 + l 



■ sm vra; 
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In this case the exact solution to the problem, is u{x, t) = e~^ * sm{iTx). We 
have performed calculations using Maple. The errors are presented in Tables 
[21 for different n- number of quadrature points (IT6|) and N- number of 
Sine-points (123|) . The table clearly exhibits an exponential decay of error 
according to the theoretical estimate. 



N 


n 


4 


8 


16 


4 


0.001074121004 






8 


0.000004530997940 


0.00000418248071 




16 


2.394152400*10-^ 


7.3086845013760*10-1° 


7.159165797001 * 10-^° 


32 


2.387505905*10-^ 


8.2307398421915*10-12 


2.609087146562 * lO-^^ 


64 




7.9836951021369*10-12 


2.917976861643 * lO'^^ 


128 






1.622297889726 * lO'^^ 


256 






1.873772451287*10-24 



Table 1: The error for x = 0.5, t 



N 


n 


32 


64 


128 
256 
512 
1024 


2.559484448336975 * lO'^^ 
2.398463850885652 * lO'^^ 
1.564339690250043*10-49 


2.3984646885635428 * 10"=^^ 

1.5716982365989911*10-49 
8.2307398421915*10-69 



Table 2: The error for x = 0.5, t 



Example 2 In this example we consider the problem as in example [I] with 
a nonlocal integral condition 



u{x, 0) 



cos s 



)u{x, s)ds 



'1 — x)x 



The results of calculation are presented in table\^for different n- number 
of quadrature points (TT6|) and N- number of Sine-points (123|) . The table 
exhibits an exponential decay of error according to the theoretical estimate. 
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n/N 


u{x,t) 


n = 4, iV = 


32 


0.5979651691 * 10"^ 


n = 8, N = 


64 


0.595184687264196427200402957709 * 10"^ 


n=16, N = 


128 


0.595184553823189342113182135931 * 10"^ 


n = 32, N = 


256 


0.595184553823189342143429937050 * 10"^ 


n = 64:, N = 


512 


0.595184553823189342143429937049 * 10"^ 


n = 128, N = 


1024 


0.595184553823189342143429937049 * 10"^ 



Table 3: Solution to the problem for x = 0.4, t = 1 
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